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(57) -M^thode pour gSn^rer un maillage hybride 
d'une formation het6rog6ne travers6e par une ou 
plusieurs discontinuites geometriques, telle que par 
exemple une formation souten^aine ou Ton afor6 un 
ou plusieurs puits. ou une formation fracturee en 
combinant des maillages structures et des mailla- 
ges non structures, dans le but de realiser des si- 
mulations en accord avec un schema numerique 
d^fini. 

Le maillage hybride dans une application a un mi- 
lieu traverse par des pults par exemple, est realise 
essentiellement en associant un premier maillage 
structure (G1) pour mailler le milieu heterogene en 
respectant ses discontinuites, des deuxiemes 
maillages structures (G2) de type radial pour mailler 



une zone autour de chaque conduit ou puits. per- 
mettant de mieux respecter des conlraintes particu- 
lieres liees aux ecoulements dans cette zone. Pour 
relier le premier maillage du milieu et les deuxiemes 
maillages de puits. on interpose entre eux des 
maillages non structures de transition (G3) et on uti- 
lise pour les realiser la technique des diagrammes 
de puissance qui est ici particulierement avanta- 
geuse en ce qu'elle permet de relier entre eux de 
maniere conforme des maillages structures non re- 
guliers. 

Applications ^ des simulateurs pour gisements 
d'hydrocarbures par exemple. 
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Description 

[0001] La pr6sente Invention concerne une m6thode pour generer un maillage hybride d'une formation heterogdne 
traversee par une ou plusleurs discontinultds g6onnetriques» dans le but de realiser par exempi© des simulations. 
[0002] La methode s'appllque particuliferement d la formation d'un maillage adapts h un gisement ou reservoir sou- 
terrain traverse par un ou plusieurs pults. ou par des fractures ou failles, dans le but d'y modeliser des deplacements 
de fluides tels que des hydrocarbures. 

Etat de la technique 

[0003] La generation de maillage est un 6l6ment cmclal pour les simulateurs de reservoir de nouvelle generation. 
Le maillage permet de ddcrire la geometrie de la structure geologique etudiee au moyen d'une representation en 
elements dlscrets dans lesquels on effectue la simulation sulvant un schema numerique approprie. Une meilleure 
comprehension des ph6nomenes physiques necessite de simuler en 3D des 6coulements polyphasiques dans des 
structures geologiques de plus en plus complexes, au voisinage de plusieurs types de singularitds tetles que les stra- 
tifications, les failles, les biseaux. les chenaux et les puits complexes. Toute cette complexite doit §tre prise en compte 
en premier lieu par le maillage qui dolt restltuer aussi fidelement que possible les informations gdologiques dans leur 
caractere h§terogene. 

[0004] La moddlisatlon de maillage a connu de grands progrfes ces demieres annees dans d'autres disciplines telles 
que i'aeronautique, la combustion dans les moteurs. la mecanlque des structures, etc. Cependant les techniques de 
maillage utilisees dans les autres domaines ne sont pas transposables telles quelles dans le monde p^troller, car les 
contraintes de m6tier ne sont pas les memes. Par exempie, en simulation de reservoir, les sch6mas numeriques sont 
construits S partir de volumes de controle afin de mieux respecter la conservation de la masse dans le cas d'equattons 
de transport de nature hyperbollque. Le maillage doit etre de type dit "block-centered" c'est-^-dire que les noeuds 
doivent etre situ^s k Tinterleur de chaque couche et les limites de chaque bloc doivent suivre I'interface entre les 
couches. Or si on ne tenait pas compte de cette contralnte. on placerait naturellement les noeuds le long des failles 
et le long des limites de stratification. Ceci aurait pour consequence de faire passer ces interfaces d travers le volume 
de contrdle utilise. La saturation, constante dans le volume de controle, ne pourrait alors respecter la discontinuity et 
les resultats seralent imprecis. II est done ndcessaire de mettre au point de nouvelles techniques mieux adaptees aux 
besoins petroliers. 

[0005] Les malllages cartSslens. couramment employes dans les simulateurs commerciaux actuals, sont mal adap- 
t6s d la resolution de ces problfemes nouveaux poses par Texploltation des gisements petroliers. Les maillages carte- 
siens, 6tant a base d'6l6ments parallelepipediques, ne permettent pas de representor de telles geometries complexes. 
[0006] On connait une methode pour generer des maillages hexaedriques structures en 3D de type dit CPG (pour 
Corner-Point-Geometry) qui respecte la geometrie des corps. On la trouve decrile dans le brevet FR 2 747 490 (US 
5 844 564) du demandeur et igalement dans la publication suivante : 

Bennis Oh. et al. « One More Step in Gocad Stratigraphic Grid G^n^ration » : Taking into Account Faults and 
Pinchouts ; SPE 35526. Stavanger, 1996. 

[0007] Ce type de maillage offre plus de souplesse que le maillage cart^sien car 11 est constitud d'dlements hexae- 
driques quelconques pouvant etre degdneres. 11 respecte rigoureusement les horizons, les failles et permet de repre- 
senter certaines inconformitSs telles que des biseaux. car sa construction s'appule sur ces elements. Toutefois ce type 
de maillage ne permet pas de resoudre toutes les complexitesgeometriques comme par exempie des maillages radlaux 
circulaires autour de pults complexes. II est possible de realiser separement le maillage du gisement et les maillages 
autour des puits mais II est difficile de representer plusieurs objets dans un meme maillage de gisement de type CPG 
a cause des problemes de raccordement Ii6s au caractere structure du maillage. 

[0008] Une autre approche est 6galement connue par laquelle on gen^re automatiquement des maillages 3D uni- 
quement^ base d'elements tetraedriques de Delaunay, avec un rafflnement radial circulaire autour des pulls. L'avan- 
tage d'une telle approche est qu'elle est entierennent automatique et ne necessite quaslment pas d'interventlon de 
rutilisateur. Toutefois cette mdthode prdsente des inconvenlents qui rendent les resultats obtenus difficilement 
exploitables : 

le nombre de mallles est en moyenne cinq fois plus grand que celul d'un maillage de type CPG pour une mdme 
structure, ce qui est tres penallsant pour les caiculs de la simulation ; 

contrairement aux maillages structures qu'on arrive alsement k visualiser. a explorer de rint6rieur et S modifier 
localement de fa^n interactive, il est trds difficile et parfols impossible de bien contrdler les malllages tetraedriques 
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a cause de leur taille et surtout du fait de leur caract^re non strublure. Ced pose probleme autant pour valider le 
maillage tfun point de vue geomdtrique que pour comprendre et valider le resultat d'une simulation sur ce type 
de maillage. 

5 [0009] D'autres approches sont egalement connues permettant de generer des maillages. notamment des maillages 
bases sur des volumes de controle gen6r6s ^ partir d'une triangulation auxquelles sont assocl^es des techniques 
d'agregatlon des triangles (ou t^tra^dres) en quadrangles permettant de redulre le nombre de mailles. Bien que des 
resultats prometteurs aient et^ obtenus avec ces nouveaux maillages. la representation precise de la complexlte geo- 
logique des reservoirs et des puits reste un sujet de recherche et de developpement. Malgre leur aspect hybride. ils 

10 restent entierement destructures et seraient done trds difflciles a gerer et a manrpuler en vrai 3D. De plus la prise en 
compte des failles reellement 3D et les pufts devies accentuerait grandement cette difficult^. 

[0010] Par la demande de brevet paralldle FR 99/15120 du demandeur. on connait une methode pour realiser en 
3D un maillage hybride sur un milieu heterogene qui comporte I'utilisation d'un premier maillage structure, de type 
CPG parexemple. et de maillages radlaux structures autourde trajectoires de puits ou conduits dessines ou importees 

15 dans le nrwdele pour mieux tenir compte des contraintes particulieres lies aux ecoulements a proximite de ces puits. 
On les combine en inserant les maillages radlaux autour du ou de chaque puits dans le maillage global du resen/oir 
apres avoir forme dans celui-ci des cavit^s assez grandes pour y former des maillages de transition non structures. 
[0011] Differentes techniques sont connues pour former des maillages non structures. On peut les realiser a base 
de poly^dres canoniques (tetraddres, pentaedres, pyramides. etc.) selon le schema numerique utilise, et I'on connaTt 

20 des solutions entierement en 3D applicables ^ ces types de maillages. 

[0012] La methode selon Tinvention vise ^ generer en 2.5D, des maillages de transition a base de polyedres de type 
VoronoV ou analogue pour y appliquer des schemas numeriques de type a volume de contrdle. que Ton ramene a un 
probleme en 2D a condition que toutes les couches des maillages de puits et du maillage du gisement pulssent se 
projeter verticalement dans un plan horizontal pour former des mailles identiques. 

25 [0013] Un tel maillage dolt respecter certaines contraintes : 

ie segment reliant les centres ou sites de deux cellules voisines doit etre orthogonal a la face commune aux deux 
cellules ; 

deux cellules voisines dolvent partager entierement la face qui les connecte. ce qui revient ^ dire que par une 
mSme face, une cellule ne peut etre connectee qu'd une seule autre cellule ; et 
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les cellules doivent dtre convexes. 
35 [0014] Ces contraintes, en 2D. imposent que : 

les aretes des polygenes frontieres ne doivent pas etre modlfiees (chacune dolt con^espondre ^ une arete d'une 
cellule creee. sans subdivision de I'arete). et que les cellules crepes doivent Stre convexes ; 

"fO - chaque cellule doit avoir un centre tel que la droite joignant les centres de deux cellules adjacentes soit perpen- 
diculaire S I'ardte partagSe par les deux cellules ; et 

les centres des cellules doivent dtre situ6s h rinterieur de leur cellule. 

45 [0015] D'autres approches entierement non structurees sont egalement connues permettant de generer des mailla- 
ges, notamment des maillages bas6s sur des volumes de controle generes a partir d'une triangulation auxquelles sont 
associ^es des techniques d'agregatlon des triangles (ou tetraedres) en quadrangles reguliers permettant de reduire 
le nombre de mailles. 

[0016] Les diagrammes de VoronoY sont formes, comme on le salt, en determinant autour d'un ensemble de sites P 
repartis dans un espace, des r6gions fomn^s d'ensembles de points qui sont plus proches de chaque site que de tout 
autre site P de I'ensemble. Cette technique de construction est decrite de fagon detailiee notamment dans les docu- 
ments sulvants : 

Aurenhammer, F., 1991 : Voronoi Diagrams : a Survey of Fundamental Geometric Data Structure. ACM Comput 
55 Sun/., 23, 346-405 

Fortune, S..1992 : Voronoi Diagrams and Delaunay Triangulations. Pages 225-265 of D.Z. Du & F.K, Hwang (eds). 
Computing in Euclidean Geometry, 2'^^edn. Lecture Notes series on Computing, vol. 4. Singapore : Worid Scien- 
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tific. 

[0017] Les cellules d'un diagramme de VoronoY respecteni les contra Intes sulvantes : Torthogonalit^ entre les centres 
et Tarete adjacente est respectee et les cellules sont convexes. II est done tout naturel de proposer de remplir la cavltS 
5 par des cellules de VoronoT» tout en cherchant a respecter les contraintes geometriques des bords. Cependant, si Ton 
cree un diagrannnne de VoronoT d partlr des points existants (extremites des aretes des polygones). (cf. Fig.6) on se 
rend compte que ce diagramnne viole les contraintes g^onnStriques, que les aretes des polygones ne font pas partie 
du diagramme et que les aretes sont coupees 

[0018] Afin d'introduire les ardtes des polygones frontieres dans le diagramme, il est possible de dupliquer les points 
10 existants. De cette fagon. le diagramme de VoronoY est oblige de passer par les aretes. Toutefois, ceci n'est pas encore 
suffisant car si les aretes sont blen dans le diagramme. elles sont encore divisees (cf . Frg.7). 

[0019] Puisque les aretes du diagramme de VoronoY correspondent aux mediatrices des segments definis par deux 
sites, on peut essayer de faire passer ces mediatrices par les extremites des aretes definies par nos polygones. Pour 
cela, pour un point extremite donne, on choisit les sites du diagramme de VoronoY le long des aretes adjacentes a ce 
IS point, a egale distance de celui-ci. Puisque les points sont d6doubtes, on est assure d'obtenir une arSte du diagramme 
de VoronoY sur Tarete de notre polygene. Dependant, si le resuttat obtenu est correct dans le cas d*un polygene rSguller 
(Fig.Sa), ce n'est pas (e cas pour un polygene non regulier (Fig. 8b). Au moment ou Ton place les demiers sites, lis ne 
sont pas forcement a la meme distance du point que les premiers que Ton avait places. 

[0020] Dans le cas g^n6ral, I'id^e d'utiliser des cellules de VoronoY pour remplir la cavity ne semble done pas ade- 
20 quate. 

La m^thode selon Tlnvention 

[0021] La methode selon Tinvention permet de simplifier la formation d'un maillage hybride adapte k un milieu he- 
25 terogene traverse par au moins une discontinuite geometrique de geometrie connue. dans le but de former un modele 
representatif d'ecoulements de fluides dans ce milieu en accord avec un schema numSrique defini. la structure du 
milieu etant connue a priori h partir de donnees disponibles acquises par des mesures in situ, par des analyses et/ou 
des interpretations d'images du milieu (des Images sismiques par exemple, dans le cas d'un gisement). 
[0022] La methode comporte la formation d'au moins un premier maillage structure (irregulier de type CPG par 
30 exemple) pour mailler au moins une partie du milieu, la formation d'au moins un deuxieme maillage structure pour 
mailler une autre partie du milieu, la formation d'au moins une cavlte entre les malllages structures englobant chaque 
discontinuity geometrique etdetaille suffisante pour que Ton puisse former au moins un maillage non structure assurant 
la transition entre les malllages structures, et la formation de chaque maillage non structure de transition en utilisant 
des diagrammes de puissance et en Imposant aux malllages de transition d'etre confomnes au maillage des parois de 
35 chaque cavite. 

[0023] La mithode s'applique par exemple ^ un milieu heterogdne ou au moins une discontinuite geometrique est 
un conduit ou puits de geometrie connue traversant le milieu, auquel cas on forme un premier maillage structure dans 
le milieu en respectant ses discontinuitds, on forme un deuxieme maillage structure de type radial autour de chaque 
puits, chaque cavity ^tant d§limit§e autour de chaque deuxieme maillage radial structure en d^sactivant des mailles 

40 de chaque premier maillage structure. 

[0024] La methode s'applique aussi par exemple a un milieu het^rog^ne oCi au moins une discontinuite geometrique 
est une fracture ou une faille traversant le milieu, auquel cas on forme un premier maillage structure et un deuxieme 
maillage structure dans des parties du milieu heterogdne, de part et d*autre de chaque fracture, en respectant leurs 
discontinuites, chaque cavite formee pour y inclure un maillage se transition non structure etant formee en desactivanl 

45 des mailles des prennier et deuxieme malllages structures, de part et d'autre de chaque faille ou fracture. 

[0025] L'utilisation qui est faite ici de la technique des diagrammes de puissance est particuli^rement avantageuse 
et commode en ce qu'elle permet de relier entre eux de manlere conforme les maillages structures non reguliers. 
[0026] L'invention concerne aussi une methode pour simuleren accord avec un schema numerique defini, revolution 
d'un processus tel que des ecoulements de fluides. dans un milieu heterogene traverse par au moins une discontinuity 

50 geometrique de geometrie connue. en accord avec un schema numenque defini. la structure du milieu etant connue 
a priori k partir de donnees disponibles acquises par des mesures in situ, par des analyses et/ou des interpretations 
d'images du milieu. Cette simulation est obtenue par la formation d'un maillage hybride constitue d'au moins un premier 
maillage structure (irregulier de type CPG par exemple) pour mailler au moins une partie du milieu, la formation d'au 
moins un deuxieme maillage structure pour mailler une autre partie du milieu, la formation d'au moins une cavite entre 

55 les malllages structures englobant chaque discontinuite geometrique et de taille suffisante pour que Ton pulsse former 
au moins un maillage non structure assurant la transition entre les maillages structures, en utilisant des diagrammes 
de puissance et en imposant aux maillages de transition d'etre confomnes au maillage des parois des dites cavites. et 
en resolvant le schema numerique dans le maillage realise sur le milieu pour modeiiser le processus. 
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[0027] Le maillage hybride global est ainsi obtenu par combinalson de plusieurs types de maillage : un malllage 
structure degisement. un maillage radial autourde chaque puits. egalement structure, etdes maillages non structures 
de transition qui assurent la liaison entre les deux types de maillage precedents. Chacun de ces maillages possede 
ses propres methodes de formation, 

[0028] En utilisant un simulateur de gisement d'un type connu, tel que par example les simulateurs ATHOS*^ ou 
SCORE™, h un gisement muni d'un maillage hybride obtenu par la mSthode, on peut realiser des simulations de 
production. 

Presentation des figures 

[002d] D'autres caractdristiques et avantages de la methode selon {'invention, apparattront a la lecture de la des- 
cription ci-aprds d'exemples non limitatifs de mise en oeuvre, en se referant aux dessins annexes oCi : 

la Fig.1 montre un exemple de maillage structure d'un resen^oir faille ; 

la Fig.2 montre un exemple de maillage radial autour d'un puits vertical ; 

la Fig.3 montre un exemple de maillage d'un puits hiorizontal ; 

20 - la Fig.4 montre en 2.5D, un exemple de reservoir maille ou Ton a menage des cavites pour des puits mailles, avant 
I'etape de creation de maillages interstitiels non structures destines a les raccorder les uns aux autres ; 

la Fig. 5 nrK>ntre un exemple de maillage en 2D avant creation du maillage de transition ; 

25 - la Fig. 6 montre un exemple de diagramme de VoronoY qui viole des contraintes geometriques ; 

la Flg.7 montre un exemple de diagramme de VoronoY dont les aretes sont divisees ; 

les Fig. 8A. 8B montrent des exemples de diagrammes de VoronoY, respective me nt dans un cas r^gulier et dans 
30 un cas quelconque avec conflit : - les Fig.9 a 14 lllustrent les definitions et propriet6s des diagrammes de puissance, 

rappeldes ci-apr^s ; 

la Fig. 15 montre un exemple de diagramme de puissance obtenu entre deux polygones ; 

35 - la Fig, 1 6 montre un exemple de maillage interstitiel en 2D forme par la methode selon invention, entre un maillage 
structure de reservoir et des maillages structures formes autour d'un puits vertical et horizontal traversant le 
reservoir : 
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la fig. 17 montre un exemple de maillage global obtenu par la methode : et 

la Fig.18 montre cinq puits pourvus chacun d'un maillage radial, integres dans un reservoir maille, par le biais de 
maillages de transition non structures a base de mailles polyedriques quelconques. 

Description deftaill6e 

[0030] La generation du maillage hybride de gisement est effectuee par etapes avec addition/soustraction de mailla- 
ges elementaires de type diff^rents. 

1) Pour rep resenter le gisennent global, on utilise un maillage structure en i, j, k d'un type connu des sp6cialistes 

50 dit CPG, comme decrit dans le brevet precite FR 2 747 490 du demandeur 

Le reservoir peut etre faille avec giissement d'un bloc par rapport a l aulre. Les horizons et les failles majeurs 
sont d'abord mod6lis6s par des surfaces continues ^ partir de donnees issues d'une interpretation d'Images sis- 
miques du sous-sol ou de donnees obtenues durant des forages (marqueurs aux puits). La structure gdologique 
est ensuite decouple en blocs failles s'appuyant sur ces surfaces. Ces blocs sont mailles individuellement; puis 

55 rassembles. Pour mailler un bloc, on commence par mailler les surfaces de bord puis on peuple I'interieur par 

interpolation transfinie des maillages surfaciques des bords. On applique alors des procedures de relaxation sur 
les surfaces de bord et h I'interieur pour harmoniser et regulariser le maillage. Le maillage ainsi obtenu respecte 
rigoureusement les horizons, les failles et permet de repr^sentercertaines inconformitds telles que des biseaux. 
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II satisfait a toutes les contraintes d'ordre geologique. 

Un tel maillage est effectue en accord par exemple avec la methode decrite dans le brevet FR 2 747 490 pr^cite. 

2] On dessine de fagon synth6tique ou on importe les trajectoires de puits. On gSn&re alors un maillage radial 
s structure autourde chaque puits pour tenlrcompte des contraintes particuli&res li^es aux 6coulements a proximite 

de ces puits. 

Dans Pexemple represented la Fig .2, le maillage structure autourd'un puits vertical est de type radial circulaire. 
C*est un nnaillage de type CPG egalement. Sa generation consiste dans un premier temps d echantltlonner un 
disque en r, 0 dans le plan horizontal. Le maillage 2D ainsi obtenu est ensuite projet^ verticalement sur tes diffe- 
10 rentes couches du maillage du reservoir. Ici, les I, J, k de la structure matricietle correspondent aux echantillonnages 

en r, 9 et z respectivement (cf. Fig. 11). 

Le maillage autour d'un puits horizontal (Fig. 3) est structure en i, j. k, du mdme type que celui du gisement a 
la difference qu'un puits ne peut pas etre faille. II est egalement obtenu en projetant verticalement sur les differentes 
couches du maillage de reservoir, un maillage 2D appartenant a un plan horizontal. 

15 

3) Purs, on insere ce maillage radial autour du ou de chaque putts dans le maillage global du reservoir. Pour cela, 
on cree dans un premier temps une cavite dans le maillage du reservoir en y desactivant toutes les mailles en 
contact avec des mailles du puits (Fig.4, 5). L'espace degage entre le maillage du gisement et le maillage du puits 
doit etre suffisant pour realiser commod6ment un maillage de transition. II peut reprdsenter par exemple environ 

20 I'equivalent de deux couches de mailles. 

4) On genere ensuite un maillage de transition non structure dans cette cavite (Fig.16« 17) pour relier le maillage 
radial structure autour du puits a celui du reservoir qui respecte au mieux ies contraintes liees au schema nume- 
rique. L'utiiisateur peut d tout moment revenird^sactiver le maillage d'un puits en reactivantles mailles de la cavite 

25 correspondante dans le maillage du reservoir. 

[0031] Pour realiser le maillage de transition on va utillser une technique connue en soi que Ton trouve decrite par 
exemple par : 

30 - Aurenhammer. F..1987 : Powerdiagrams : properties, algorithms and applications. S/AM J. Comput., 16(1). 78-96. 
dont on va rappeter ci-aprds quelques notions pour la clarte de Texpose. 
Point pondere 

35 

[0032] Soit p un point desi^ auquel on associe un poids w. On appelle point pondere et on note p le couple (p. vi^ 
ainsi forme. Un point pondere peut aussi etre interprets comme un cercle C de centre p et de rayon Jw. 

Notion de puissance 

40 

[0033] On designe par Puis^p^^ = PuiSfs^x) = Ip^I^ - puissance d'un point x par rapport au point pondere p ou par 
rapport au cercle C. |px|2 represente ici la distance euclidienne entre p et x. Geometriquement. si x est a I'exterieurdu 
cercle et si t est le point d'intersection d'une tangente a C passant par x, Puis^^^ = |xf|2 (cf. Fig. 9). 

45 Puissance d'un point par rapport a un cercle 

[0034] Soit C le cercle de centre c et de rayon r, d'equation : 

so 0{x) = 0 ou C(x) = xc^ - r 

C(x) est la puissance du point x par rapport ^ C. 
Diagramme de puissance 

55 

[0035] Soit S un ensemble de points ponderSs. A chaque point p de S, on associe une region /?(p), potentlellement 
vide, definie par : 
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R{p)^{xsSx'/ P«^p.,o < Piiis^^ ,,,, P. e 5 - {p}} 

[0036] Get ensemble de regions, ou encore cellules est appele diagramnne de puissance de S. 
Axe radical de deux cercles 

[0037] L*axe radical de deux cercles est la droite formee par les points ayanl la meme puissance vis d vis des deux 
cercles. En particulier si les deux cercles ont une intersection non nulle I'axe radical est la droite passant par Tinter- 
section (Flg.10). 

Utilisation des diagrammes de puissance 

[0038] Les diagrammes de puissance se pretent bien a la resolution du probleme pose car lis ont I'avantage de 
possSder les proprietes suivantes : les cellules qui te composent sent convexes, chaque cellule poss^de un centre 
(site) et la droite reliant les centres de deux cellules adjacentes est perpendtculaire S Tarete commune aux deux cellules, 
proprietes qui satisfont done a une partie des conditions posees. 

[0039] On souhaite en outre que les cellules du diagramme de puissance respectent les aretes des polygones fron- 
tieres. II faut done que les aretes des polygones soient des aretes du diagramme de puissance. Pour r^aliser un tel 
diagramme, on impose comme contrainte que les aretes des polygones soitles arStes d une triangulation de Delaunay 
couramment utilis^e dans le domaine des maillages. Cette contrainte estfacilement realisable car on determine les 
polygones frontiSres lors de la creation de la cavite. Une fois cette condition etablie, il est possible comme on va le 
voir, de creer un diagramme de puissance respectant les contraintes. 

Diagramme de puissance contraint : respect des aretes donnees 

[0040] Soit (S) un ensemble de n aretes dans le plan appartenant a une triangulation de Delaunay {T^) (cf. Fig. 11). 
[0041] Soit (S*) I'ensemble des aretes duales de (S)dans te diagramme de VoronoY associ6 d (T^)(cf, Fig. 12). 
[0042] Soit a, et bf, i = 1..n, les extremites de I'arete S/. On considers son arSte duale S*/ dans le diagramme de 
Voronoi associe a (TJ). Elle definit un ensemble de cercles (C;) de centres Cf pris sur 8% et de rayon || c/aj || = || Cjbj ||. 
Tout cercle C/ de {Cj) porte done a/ et 6/ et. est un cercle vide (Fig. 13). On considere C\ et deux cercles de {Cj). 
D'apr6s ce qui precede, a, et />,- sont sur Cy et C^j, et bf sont done les points d'intersection de C\. et C^j et se situent 
sur leur axe radical. En termes de puissance, on peut done ecrire : 

C]{a,} = Cf(a,) = cUbf) = Cf(b,) = 0 

[0043] Par construction, quelque soit Cy de {Cpj = L./i^yV /, I'interieur de Cy est vide et ne contlent done ni a;, ni 6,-. 
On peut encore ecrire : 

y/. C/a/)^0 

vy. Cj(b,)^o 

[0044] De 1^, Bi appartlent a une arete du diagramme de puissance, portee par Taxe radical de C] et de , bj 

appartient a une arete du diagramme de puissance, portee par I'axe radical de C/ et de C^- (Fig. 14). 

[0045] Comme par definition, les cellules d'un diagramme de puissance sont convexes, on en d6duit que le segment 

L Bf b^ = S/fait partie du diagramme de puissance. 

[0046] Ceci est vrai pour tout S, de (S). On peut done construire un diagramme de puissance contraint contenant 
I'ensemble des segments de (S). Le degre de liberie pour construire un tel diagramme est grand car la seule contrainte 
est de choisir deux centres par arete duale S*,- pour tout de (S). 

[0047] On remarque que, si S*f est r6duite a un point. I'ensemble {Cf) est redult ^ un cercle de centre C; = S*/. II faut 
done utiliser dans un tel cas un centre virtuel. infinlment proche du premier pour ddfinir un second cercle. 
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Construction du maillage de transition ^ i'aide des diagrammes de puissance 

[0048] Pourcreerun diagrannme de puissance, maillage de transition, on commence parconstruire une triangulation 
de Delaunay des extremitds sommets des arStes des polygenes fronti^res. On verifie que toutes les aretes initiates 

s sent blen des aretes de la triangulation obtenue. Ceci fait, on construit le diagramme de VoronoY associS aux sommets 
des polygones. Ce diagramme de VoronoY d^finit pour chaque arete A des polygones. une ardte duale A ' 
[0049] Par la suite, pour chaque arete A, on positionne deux sites du diagramme de puissance final sur A Puisque 
les sites doivent etre des points ponderes. on associe au site x situe sur A' un poids equivalent d la distance \\ax\\^ oCi 
a est une extremite de A. Le diagramme de puissance obtenu a partir de ces sites respecte alors les aretes des 

10 polygones. La Fig. 16 est un exemple de diagramme de puissance obtenu entre deux polygenes avec un seul site par 
arete. 

[0050] Une fois ce diagramme obtenu, toutes les contraintes de depart sont maintenant respectees sauf a bien 
verifier que chaque site est h t'interieur de sa cellule. Par construction, on a positionne les sites sur les aretes duales 
sans en pr^clser la mani^re. On a jusqu'a present montr6 qu'il existait un diagramme de puissance con-ect d^s lors 
IS que les sites dtaient choisis sur A *, mais les possibilitds dans le choix de ces sites sont multiples. II est possible 
toutefois de garantir que si touts arete A intersecte son arete duale A\ il existe au moins un diagramme de puissance 
respectant les polygones fronti^res et dont les sites sont ^ I'interieur de leur cellule. 

R^sultats 

20 

[0051] La Fig. 16 montre deux exemples de maillage 2D. Le premier exemple correspond a une situation de puits 
vertical quant ^ la seconde elle represente le cas d'un puits horizontal. On peut verifier sur ces exemples que les 
cellules g6nerees pour jolndre les maillages des puits et du gisement respectent les contraintes enoncees dans le 
problems : les cellules sont convexes, les sites sont ^ I'interieur de leur cellule, les aretes des polygones (bords des 
2S maillages) sont respectees. et la droits joignant les sites de deux cellules adjacentes est perpendiculaire a t'arete 
commune aux deux cellules. 

[0052] L*6tape 2D permet dans un premier temps de generer les maillages de transition et ainsi de connecter les 
maillages de puits avec le maillage de gisement. Pour obtenir des maillages en volume, on transforme les maillages 
initiaux en maillages 2D, et on applique aux resuttats en 2D un traitement inverse, permettant d'obtenir un maillage 
30 global final en 2.5D. 

[0053] On a consid6r6 le cas ou les deuxiemes maillages structures etaient de type radial et formes autour de puits 
traversant le milieu, avec delimitation de cavites autour de chaque deuxieme maillage pour y inclure un maillage de 
transition. II est clair cependant que la methods s'applique au maillage d'un milieu presentant d'autres types de dis- 
continuity g6om6trique, te! que par exemple un gisement soulerrain traverse par des fissures. Dans un tel cas, on 
35 forme des maillages structures dans des parties du milieu de part et d'autre des fissures et Ton delimite des cavites 
en desactlvant des mallles des deux maillages structure autour de chaque fissure, assez grandes pour y former des 
maillages de transition non structures en utilisant des diagrammes de puissance et en imposant aux maillages de 
transition d'etre conformes au maillage des parois de chaque cavite. 

40 

Revendications 

1, Methode pour generer un maillage hybride adapte d un milieu heterogene tel que le sous-sol, traverse par au 
moins une discontinuite geometrique de geometrie connue, dans le but de former un modele representatif d'ecou- 

4S lements de fluides dans ce milieu en accord avec un schema numerique deflni. la structure du milieu etant connue 

a priori k partir de donnees disponibles acquises par des mesures in situ, par des analyses et/ou des interpretations 
d'images du milieu, comportant la formation d*au moins un premier maillage structure pour mailler au moins une 
partle du milieu, la formation d'au moins un deuxieme maillage structure pour mailler une autre partie du milieu 
contenant la discontinuite geometrique. la formation d'une cavite entre le premier maillage structure et chaque 

so deuxieme maillage structure englobant chaque discontinuite geometrique et de taille suffisante pour que Ton pulsse 

former au moins un maillage non structure assurant la transition entre les maillages structures, caracterisee en 
ce que Ton forme chaque maillage non structure de transition en utilisant des diagrammes de puissance et en 
imposant aux maillages de transition d'etre conformes au maillage des frontieres de chaque cavite. 

55 2. Methode selon la revendication 1, appliquee e un milieu heterogene oij au moins une discontinuite geometrique 
est un conduit ou puits de geometrie connue traversant le milieu, caracterisee en ce que chaque deuxieme maillage 
structure est de type radial et chaque cavite est deiimitee autour de chaque deuxieme maillage radial structure en 
desactlvant des mallles de chaque premier maillage structure. 
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Methode selon la revendication 1 ou 2, appliquSe h un milieu heterog^ne ou au molns une discontinuite gSome- 
trique est une fracture ou une faille traversant le milieu heterogene. caracterisee en ce que Con fomie un premier 
maillage structure et un deuxi^me maillage structure dans des parties du milieu heterogene. de part et d'autre de 
chaque fracture, en respeclant leurs discontinuit^s. chaque cavite formee pour y inclure un maillage se transition 
non structure etant formee en desactivant des mailles des premier et deuxieme maillages structures, de part et 
d'autre de chaque fracture. 

M6thode selon I'une des revendicatlons precedentes. caracterisee en ce que Ton impose aux aretes polygonales 
formant les fronti^res de chaque cavitS. d'Stre les ardtes d'une triangulatlon de type Delaunay. 

Methode pour simuler en accord avec un schema num6rique d^finl, revolution d*un processus tel que des ecou> 
lements de fluides dans un milieu heterogdne traverse par au motns une discontinuite geometrique de geometrie 
connue. la structure du milieu etant connue a priori k partir de donnSes disponibles acquises par des mesures in 
situ, par des analyses et/ou des Interpretations d'images du milieu, comportant la fonnation d'un maijiage hybride 
constitue d'au molns un premier maillage structure pour mailler au moins une partie du milieu, la formation d'au 
moins un deuxieme maillage structure pour mailler une autre partie du milieu, la formation d'une cavite entre le 
premier maillage structure et chaque deuxieme maillage structure englobant chaque discontinuite geometrique et 
de tailte suffisante pour que Ton puisse former au moins un maillage non structure assurant la transition entre les 
maillages structures, caracterisde en ce que : 

Ton construit les maillages non structures de transition en utilisant des diagrammes de puissance et en impo- 
sant aux maillages de transition d'etre conformes au maillage des parois des dites cavltes ; et 

on r6sout le schema num^rique dans le maillage realise sur le milieu pour modellser le processus. 

Methode selon Tune des revendicatlons prec6dentes, caracterisee en ce que chaque premier nnaillage structure 
est un maillage non r^gulier, de type CPG par exemple. 
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